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Abstract 


Laser ablation (LA) and spark discharge (SD) techniques are commonly used for 
nanoparticle (NP) formation. The produced NPs have found numerous applications 
in such areas as electronics, biomedicine, textile production, etc. Previous studies pro- 
vide us information about the amount of NPs, their size distribution, and possible ap- 
plications. On one hand, the main advantage of the LA method is in the possibilities 
of changing laser parameters and background conditions and to ablate materials with 
complicated stoichiometry. On the other hand, the major advantage of the SD techni- 
que is in the possibility of using several facilities in parallel to increase the yield of 
nanoparticles. To optimize these processes, we consider different stages involved and 
analyze the resulting plasma and nanoparticle (NP) parameters. Based on the per- 
formed calculations, we analyze nanoparticle properties, such as mean size and mean 
density. The performed analysis (shows how the experimental conditions are connect- 
ed with the resulted nanoparticle characteristics in agreement with several previous 
experiments. Cylindrical plasma column expansion and return are shown to govern 
primary nanoparticle formation in spark discharge, whereas hemispherical shock de- 
scribes quite well this process for nanosecond laser ablation at atmospheric pressure. 
In addition, spark discharge leads to the oscillations in plasma properties, whereas 
monotonous behavior is characteristic for nanosecond laser ablation. Despite the dif- 
ference in plasma density and time evolutions calculated for both phenomena, after 
well-defined delays, similar critical nuclei have been shown to be formed by both 
techniques. This result is attributed to the fact that whereas larger evaporation rate is 
typical for nanosecond laser ablation, a mixture of vapor and background gas deter- 
mines the supersaturation in the case of spark. 


Keywords: Nanoparticles, laser ablation, plasma, spark discharge, synthesis, model- 
ing, size distribution, nucleation, aggregation 
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1. Introduction 


Modern nanotechnology includes several promising areas such as nano-optics, nano-photon- 
ics, nanochemistry, nanobiology, and nanomedicine. During the past decade, we have 
witnessed a tremendous growth of nanoparticle applications that require particles of different 
materials with different size, dispersion, shape, and morphology. As a result, the development 
of new nanoparticle synthesis methods is particularly important. Among the rapidly im- 
merged techniques, plasma-based synthesis has a number of advantages being both rather 
simple and allowing unique and well-controlled formation of nanoparticles. 


In particular, plasmas created by both pulsed laser ablation and by spark discharges can be 
used for nanoparticle synthesis. That is why these two techniques have attracted particular 
attention and resulted in numerous experimental and theoretical investigations. On one hand, 
the main advantage of laser ablation method, as was demonstrated in these studies, is the 
possibility to preserve target stoichiometry. On the other hand, spark discharge allows one to 
produce a very large amount of nanoparticles by using parallel multidischarge set-up. 
However, the main physical mechanisms involved in these processes stay partly puzzling. 
That is why additional studies based on a detailed comparison of both methods are required 
for the determination of their main similarities and differences. 


Starting from the early 1990s, laser ablation (LA) has been intensively studied first for long 
laser pulses and then for much shorter ones [1-2]. A number of experiments revealed that laser 
interactions with solid targets lead to the formation of nanoparticles. Furthermore, if femto- 
second laser is used, an explosive ejection of a mixture of clusters and atoms was both 
theoretically predicted and experimentally observed [3]. It was found that the produced 
nanoparticles demonstrated either plasmonic or photoluminescent properties, as well as a 
capacity of field amplification. These properties are particularly interesting biomedical 
applications, such as imaging and photodynamic therapy. It should be noted that absence of 
incompatibility with biological tissues is crucial for further development of most of the 
applications where nanoparticles are used in vivo. However, toxicity is hard to avoid in 
traditional chemical methods. In addition, the stability of these nanoparticles is still not high 
enough. It was demonstrated, fortunately, that nanoparticles produced by laser ablation are 
better suitable for biomedical applications, in particular, when they are produced in liquids [4]. 
This advantage made LA a unique tool for nanoparticle synthesis. 


In order to elucidate the physical mechanisms of LA, many analytical and numerical investi- 
gations were performed [5 -11]. In vacuum, self-similar adiabatic models with condensation 
were proposed [10-11]. In the presence of a gas, only either very low or high background 
pressure was used in most of the models for simplicity. Shock waves were shown to be 
produced during the plume expansion in a high-pressured background gas [3]. In this case, a 
system of Navier-Stocks equations well describes the first 1-2 us of the plasma plume expan- 
sion. It should be noted that such models are invalid at the later stages. To solve this issue, 
hydrodynamic calculations were switched to the direct Monte Carlo simulations where no 
such hypothesis is used at 1-2 us after the laser pulse [11, 12]. Recently, such approaches as 
molecular dynamics (MD), hydrodynamics (HD), and combinations with the direct simulation 
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Monte Carlo method (DSMC) were proposed for picosecond and femtosecond laser interac- 
tions [13-14]. The initial stage of laser ablation process can be examined by using either 
hydrodynamic models or atomistic simulation. On one hand, the main advantage of numerical 
hydrodynamics is in the calculation rapidity and in the possibility to reach rather larger scales 
[10, 13]. Atomistic simulations, on the other hand, are not based on equilibrium assumptions 
and can more easily provide size distributions of nanoparticles. In particular, two-temperature 
molecular dynamics simulations (TTM-MD) were performed for femtosecond laser ablation 
of metals [14, 15]. 


In the presence of a sufficiently high-pressured background environment, such as atmospheric 
pressured gas or a liquid, diffusion-driven nucleation and aggregation processes start playing 
an important role (Figure 1) at longer delays [16]. 
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Figure 1. Nanoparticle formation and growth schematics. 


However, despite a growing interest in LA, some of its basic mechanisms remain not enough 
understood. The challenge is that ablation processes strongly depend on the ensemble of such 
parameters as temporal pulse width and shape, on laser wavelength, on the size of the laser 
spot, laser intensity, repetition rate, as well on the target material and background conditions. 
[17] 


In spark discharges (SD), a typical set-up consists of two electrodes connected to a charged 
capacitance [16]. When high enough voltage V, is applied, a so-called streamer is formed first. 
Once it reaches the opposite electrode, plasma breakdown takes place, followed by the 
streamer transition to an expanding plasma column. In this column, Joule heating of both 
plasma and electrodes take place. In addition, electrodes are bombarded by energetic ions that 
induce sputtering, which can be considered to be similar to laser ablation. If background gas 
is present, rapid thermalization of the sputtered material leads to primary nanoparticle 
formation that can then grow by collisions and form larger particles or aggregates. 


Numerical modeling of spark discharge consists of several steps with rather different time 
scales [16]: (i) streamer formation and propagation between electrodes; (ii) streamer-to-spark 
transition; (iii) gas heating and cylindrical expansion, (iv) electrode evaporation and erosion; 
(v) nanoparticles formation. Nanoparticle formation, here as in the case of laser ablation in the 
presence of a gas or a liquid, includes nucleation and collisional growth (Figure 1). 
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This work is aimed at the better understanding of the mechanisms involved in nanoparticle 
formation by laser ablation and by spark discharge. First, laser ablation is considered by using 
both atomistic and hydrodynamic numerical methods. In particular, mechanisms of nanopar- 
ticle formation and the corresponding conditions are analyzed. Then, attention is focused on 
the role of the background environment and its role in nanoparticle formation. These results 
are used to explain several recent experimental results. 


Second, spark discharge is investigated. Plasma properties and conditions required for 
nanoparticle formation are examined. Finally, we compare laser ablation and spark discharge 
as promising methods of nanoparticle formation. 


2. Laser ablation 


To examine ultra-short, laser-ablated plume dynamics and nanoparticle evolution under 
realistic experimental conditions and to account for the fact that the ablated plume contains 
several components, DSMC calculations of the plume dynamics are first performed in the 
presence of an inert background gas (Ar) with pressure P = 300 Pa. The initial conditions are 
set based on the parameterization of the MD results obtained at a delay of 200 ps after the 
beginning of the laser pulse (100 fs, 800 nm) [16]. 


Figure 2 shows separately the density of atoms and clusters for two different delays after the 
laser pulse. Here, larger clusters were initially at the back of the plume. The obtained results 
demonstrate that plume front starts experiencing a pronounced deceleration and practically 
stops at both plume- and gas-dependent delay (here, ~10us). Theoretically, the initial expan- 
sion stage is described by a so-called blast-wave (or, shock-wave) model when shock waves 
are degenerated. 


The corresponding nanoparticle size distributions are presented in Figure 3. One can see that 
after a sufficient delay, a peaked distribution appears instead of a decreasing function. This 
effect can be explained by collisional growth that is described by the general rate equation 
having typically log-normal solutions. The amount of sufficiently large nanoparticles formed 
at such short delays is rather small and cannot explain the finally observed size distributions. 


Longer stage includes plume mixing with the background followed by the rapid thermaliza- 
tion and a much more enhanced particle formation. Then, plume species are thermalized and 
a diffusion-driven regime enters into play. 


3. Spark discharge 


In our model developed for spark discharge (SD), the following stages are considered: (i) 
streamer formation and propagation between electrodes; (ii) streamer-to-spark transition; (iii) 
gas heating and cylindrical expansion; (iv) electrode evaporation and erosion; (v) nanoparticles 
formation [16]. 
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Figure 2. Calculated plume dynamics for Ni expansion in Ar gas at 300 Pa, (a) density snapshot for atoms att = 0.55 us, 
(b) the same for clusters at t = 0.55 us; (c) the same for atoms at t = 10 us, (d) the same for clusters at t = 10 us. 


Streamer formation is described numerically by using a system of drift-diffusion equations 
together with Poisson equation for electric potential with a particularly chosen set of boundary 
conditions. When streamer reaches the opposite electrode, electron emission increases 
dramatically, so that the streamer is transformed in a conductive plasma column. The oscilla- 
tions of the electric charge Q in the corresponding circuit with the total resistance Rs is 
described by the Kirchhoff’s voltage law. 


The oscillating behavior of the discharge is presented in Figure 4, defining the properties of 
the following plasma column. According to this solution, both electrodes play a role at different 
delays leading to both evaporation and erosion of the electrodes. When polarity switches, a 
crater is formed at the surface of one of the electrodes due to both evaporation and erosion. 


The erosion flux jak is formed due to two main processes: (i) thermal evaporation caused by 
Joule heating, and (ii) sputtering due to ion bombardment. Similarly to nanosecond laser 
ablation, evaporation flux ae is described by the Hertz-Knudsen equation, where surface 


temperature is calculated as follows [16] 
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Figure 3. Size distributions calculated by using MD-DSMC model in the presence of 300 Pa of Ar at different time de- 
lays. 
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Figure 4. Typical time evolutions of voltage and electric current during a “single” spark event obtained in N, for C =8 
nF; L = 0.77 uH; Ry =1Q. 
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where s corresponds to the cathode material (solid), c is the specific heat of cathode material, 
pis the density, T, is the temperature, o, is the electric conductivity, x, is the thermal diffusivity 
coefficient, k is the Boltzmann constant, a is the sticking coefficient for vapor molecules onto 
the surface, P, is the equilibrium vapor pressure, and P,,,,,is the hydrostatic pressure of gas 
applied on the surface. 


The sputtering yield Y Y and flux jsut is given by 


4 


L Sd nf 4mM E 
An? (m + MÝ U urf ‘ (2) 


‘sput ‘sput 


. F i 
I surf = LY J surf = J surf g I surf 7 


where @,,,,;is a factor function of m/M, m and M are the atomic weights of cathode material and 
incident particles, respectively, E, is the bombarding energy, U.,,;is the surface-binding energy, 
and j, is the flux density of bombarding ions. The solution of Eqs (1-2) yields cathode erosion 
flow as a function of time. At the same time, plasma column also gains energy by Joule heating. 


Once the amount of the ejected material is calculated, plasma dynamics is modeled by using 
Navier-Stokes equations [17]. The corresponding equations contain, in particular, Joule 
heating term, which determines plasma heating occurring mostly near its axis, where plasma 
pressure initially arises to several atmospheres. Figure 5 shows typical axial temperature 
evolution. After a delay of 0.5 us, which corresponds to plasma expansion, pressure drops 
back to the values around the atmospheric pressure. Plasma temperature remains high during 
all the discharging process and drops below electrode boiling point only after 0.1-0.5 ms. 
During cooling, gas density increases. Lager density and smaller temperature provide 
conditions required for nanoparticle formation. 
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Figure 5. Typical time evolutions of the plasma temperature obtained for Ar for gap of 1 mm, R =1 Q, C=8nF, and L= 
0.77 WH. 
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4, Spark discharge vs nanosecond laser ablation 


Typically, laser energy absorption leads to the target heating and thermal evaporation in the 
nanosecond laser ablation of metals. Plasma expansion stage is much longer than the evapo- 
ration stage, on the order of ~10 laser pulse temporal widths (~300 ns). Here, plasma electrons 
gain energy from laser radiation by inverse Bremsstrahlung effect, so that ionization takes 
place. Then, a so-called blast wave model describes hemispherical expansion as follows [16]: 
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where R, is the shock wave position, € is a constant depending on the plume specific heat 
capacity, Eois the initial internal energy of vapor plume, Tọ is the initial temperature, and y is 
the adiabatic coefficient. This model is valid, if the external shock wave is present and if the 
mass of the background gas surrounding the shock wave is larger than the mass of the ablated 
material. During the expansion, both plume temperature and density rapidly decay. As a 
result, the condition of supersaturation is realized at a certain delay leading to the nanoparticle 
formation. 


Time evolutions of plasma properties are obtained in [16]. For spark discharge, when gas 
density decreases, plasma resistivity drops down, thus preventing further energy absorption 
from the electric current. During the expansion, plasma cools down and gas goes back to the 
axis. These processes are repeated when the energy absorption is high enough for efficient 
expansion. 


In the case of nanosecond laser ablation, laser energy absorption by the plume takes place. As 
aresult, the ejected plasma is hot, but the energy input is limited in time by laser pulse duration. 
After the end of the laser pulse, the ablated plume expands and cools down. The critical nuclei 
sizes are calculated as a function of time (Figure 6). 


The obtained results clearly demonstrate that, despite rather different plasma sources, the 
behavior of particle formation in spark discharge appears to be very similar to the one in laser 
ablation. The difference in time delays corresponds to the transition of the system to the 
supersaturated state. Primary NPs formed by nucleation can then evolve due to evaporation, 
condensation, and/or growth. 


Nanoparticle Formation by Laser Ablation and by Spark Discharges — Properties, Mechanisms, and... 
http://dx.doi.org/10.5772/61303 


30 


— SD 
25 —LA 


Critical radius [nm] 


0 r T r T r T r T r T r 1 
0.0 0.5 1.0 1.5 2.0 2.5 3.0 


Time [ys] 


Figure 6. Typical time evolution of primary critical clusters size formed by spark discharge and by laser ablation. 


Diffusion-driven nucleation leads to the formation of nucleus, whose size is controlled by the 
free energy as follows [18]: 


AG(n,c) =—nkTIn(c/c,, ) + 4Aran?o, (6) 


where k is the Boltzmann constant; T is the temperature in Kelvins; a is the effective radius; ca 
is the equilibrium concentration of atoms/monomers; n is the number of atoms/monomers in 
the nuclei; and o is the effective surface tension. The peak of the nucleation barrier corresponds 


to the critical cluster size 
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The production rate of the supercritical nuclei is given by 


v(t)=K e ap SO] (8) 


As a result, narrow size distributions are produced. The formed particles can collide and 
aggregate. These collisional processes are described by simplified Smoluchowski equations 
[18]. Laser-induced fragmentation is neglected here. 


Several calculations are performed to study time-evolution of the size distribution for multiple 
pulse cases. Because saturation is rather high, monomer radius here is as small as a = 1.597 
and the radius of critical nuclei is recalculated at all the time-steps. 
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Figure 7. Calculated size distribution obtained for 10, 100, 1000, and 2000 pulses. Here, laser frequency is 1 kHz, gold 
solution in water is considered with a = 1.59 m. 


The obtained results (Figure 7) show that when several pulses are applied, these small nuclei 
grow for pulse number up to 10° because tiny particles grow easier than larger ones. Here, 
collisions with atoms dominate in the growth process. The obtained results show that further 
increase in the number of pulses affects particle size distribution only slightly. This effect takes 
place if laser frequency is not too high, so that the ablated material has time to diffuse and 


concentration does not grow near the target. As a result, the mean radius can remain rather 
small (in the nanometer range, smaller than ~3 nm here). 


The criterion of the catastrophic nucleation due to thermalization can be derived based on the 
inequality n, < 1. This means that 


N, =| 8700 /3k,Tin(c/c,,)] ~1. 9) 


Typically, at the beginning of the diffusional expansion stage, the condition 
In(S)=In(c / Coq) = 870070 /3kT 21 is valid, so that this mechanism prevails in the nanoparticle 
formation. Finally, at the last stage, ripening or sintering of the created particles eventually 
inters into play if the background density is sufficiently high. For instance, this process occurs 
in liquids, in particular, in the absence of the surface passivation by additional chemicals. 


5. Conclusions 


In this chapter, recent advancements in the modeling of laser ablation and spark discharges 
are summarized. In addition, we have compared processes leading to nanoparticle formation 


in laser ablation and in spark discharges. In particular, a comparison of the influence of plasma 
properties on nanoparticle formation has been performed. 
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First, mechanisms of nanoparticle generation have been investigated for femtosecond laser 
interactions in the presence of a background environment. The obtained calculation results 
have demonstrated the long time-evolution of plume species involves nucleation and growth 
determining the final size distribution that tends to be a limited one with the increase in laser 
intensity. The obtained results explain several experimental observations including both 
longer time dynamics of nanoparticles and size distributions. Then, conditions are formulated 
for catastrophic nucleation to become the main mechanism of nanoparticle formation as a 
result of thermalization and collisions among the species in the presence of a background 
environment. The produced nanoparticles can be collected, form a colloid, or can be deposited 
at a substrate forming nanostructures. Therefore, the presented study is of interest for many 
applications where both metallic nanoparticles and nanostructures are used in nanophotonics, 
plasmonics, medicine, biological sensing, textile industry, and other promising fields. 


Furthermore, cylindrical plasma column expansion has been shown to govern primary 
nanoparticle formation in spark discharge, whereas hemispherical shock describes quite well 
this process for nanosecond laser ablation at atmospheric pressure. In addition, spark dis- 
charge leads to oscillations in plasma properties, whereas monotonous behavior is a charac- 
teristic for nanosecond laser ablation. 


Despite the difference in plasma density and time evolutions calculated for both phenomena, 
after well-defined delays, similar critical nuclei have been shown to be formed by both 
techniques. This result can be attributed to the fact that whereas larger evaporation rate is 
typical for nanosecond laser ablation, a mixture of vapor and background gas determines the 
supersaturation in the case of spark. 
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